

*****************************************************
* (1.0) PPR Robustness Checks
*****************************************************
	global controls1 	"m_year_birth i.educ age_marriage"
	global controls2 	"grain agric_production gdp pop_total_agric"
	global controls3 	"U5MR_5yr_ave"
	
	************************************************
	* (1.1) Robustness Regressions
	************************************************
	mat robustness = results 
	
	forval spec = 1/7 { 
	
		********************************************
		* No Controls
		********************************************
		else if `spec' == 1 {
			xi: reg delivered years_since_LLFpre8-years_since_LLFpre2 years_since_LLFpost0-years_since_LLFpost7 ///
				pre_8_son-pre_2_son post_0_son-post_7_son ///
				p1_LLFpre8-p1_LLFpre2 p1_LLFpost0-p1_LLFpost7  ///
				p2_LLFpre8-p2_LLFpre2 p2_LLFpost0-p2_LLFpost7 ///
				p2_LLFpre8_son-p2_LLFpre2_son p2_LLFpost0_son-p2_LLFpost7_son ///
				p3_LLFpre8-p3_LLFpre2 p3_LLFpost0-p3_LLFpost7 ///
				p3_LLFpre8_son-p3_LLFpre2_son p3_LLFpost0_son-p3_LLFpost7_son ///
				p_1 p_2 p_3 hasson  p2_son p3_son i.FE_t i.FE_prov ///
				if year!=1967 & year !=1968 & years_since_LLF<8 & rural==1, cluster(prov)
		}
		
		************************************************
		* Individual Characteristic Controls Only
		************************************************
		else if `spec' == 2 {
			xi: reg delivered years_since_LLFpre8-years_since_LLFpre2 years_since_LLFpost0-years_since_LLFpost7 ///
				pre_8_son-pre_2_son post_0_son-post_7_son ///
				p1_LLFpre8-p1_LLFpre2 p1_LLFpost0-p1_LLFpost7  ///
				p2_LLFpre8-p2_LLFpre2 p2_LLFpost0-p2_LLFpost7 ///
				p2_LLFpre8_son-p2_LLFpre2_son p2_LLFpost0_son-p2_LLFpost7_son ///
				p3_LLFpre8-p3_LLFpre2 p3_LLFpost0-p3_LLFpost7 ///
				p3_LLFpre8_son-p3_LLFpre2_son p3_LLFpost0_son-p3_LLFpost7_son ///
				p_1 p_2 p_3 hasson  p2_son p3_son $controls1 i.FE_t i.FE_prov ///
				if year!=1967 & year !=1968 & years_since_LLF<8 & rural==1, cluster(prov)
		}
		
		********************************************
		* Individual and Provincial Controls
		********************************************
		else if `spec' == 3 {
			xi: reg delivered years_since_LLFpre8-years_since_LLFpre2 years_since_LLFpost0-years_since_LLFpost7 ///
				pre_8_son-pre_2_son post_0_son-post_7_son ///
				p1_LLFpre8-p1_LLFpre2 p1_LLFpost0-p1_LLFpost7  ///
				p2_LLFpre8-p2_LLFpre2 p2_LLFpost0-p2_LLFpost7 ///
				p2_LLFpre8_son-p2_LLFpre2_son p2_LLFpost0_son-p2_LLFpost7_son ///
				p3_LLFpre8-p3_LLFpre2 p3_LLFpost0-p3_LLFpost7 ///
				p3_LLFpre8_son-p3_LLFpre2_son p3_LLFpost0_son-p3_LLFpost7_son ///
				p_1 p_2 p_3 hasson  p2_son p3_son $controls1 $controls2 i.FE_t i.FE_prov ///
				if year!=1967 & year !=1968 & years_since_LLF<8 & rural==1, cluster(prov)
		}
		
		*****************************************************
		* Individual + Provincial + Mortality Controls
		*****************************************************
		else if `spec' == 4 {
			xi: reg delivered years_since_LLFpre8-years_since_LLFpre2 years_since_LLFpost0-years_since_LLFpost7 ///
				pre_8_son-pre_2_son post_0_son-post_7_son ///
				p1_LLFpre8-p1_LLFpre2 p1_LLFpost0-p1_LLFpost7  ///
				p2_LLFpre8-p2_LLFpre2 p2_LLFpost0-p2_LLFpost7 ///
				p2_LLFpre8_son-p2_LLFpre2_son p2_LLFpost0_son-p2_LLFpost7_son ///
				p3_LLFpre8-p3_LLFpre2 p3_LLFpost0-p3_LLFpost7 ///
				p3_LLFpre8_son-p3_LLFpre2_son p3_LLFpost0_son-p3_LLFpost7_son ///
				p_1 p_2 p_3 hasson  p2_son p3_son $controls1 $controls2 $controls3 i.FE_t i.FE_prov ///
				if year!=1967 & year !=1968 & years_since_LLF<8 & rural==1, cluster(prov)
		}
		
		*****************************************************
		* Prov*Trends
		*****************************************************
		else if `spec' == 5 {
	
			xi: reg delivered years_since_LLFpre8-years_since_LLFpre2 years_since_LLFpost0-years_since_LLFpost7 ///
				pre_8_son-pre_2_son post_0_son-post_7_son ///
				p1_LLFpre8-p1_LLFpre2 p1_LLFpost0-p1_LLFpost7  ///
				p2_LLFpre8-p2_LLFpre2 p2_LLFpost0-p2_LLFpost7 ///
				p2_LLFpre8_son-p2_LLFpre2_son p2_LLFpost0_son-p2_LLFpost7_son ///
				p3_LLFpre8-p3_LLFpre2 p3_LLFpost0-p3_LLFpost7 ///
				p3_LLFpre8_son-p3_LLFpre2_son p3_LLFpost0_son-p3_LLFpost7_son ///
				p_1 p_2 p_3 hasson  p2_son p3_son culturalrev $controls1 $controls2 $controls3  i.year prov##c.year i.FE_prov ///
				if year!=1967 & year !=1968 & years_since_LLF<8 & rural==1, cluster(prov)
		}	
	
		*****************************************************
		* Including CR years
		*****************************************************
		else if `spec' == 6 {
			xi: reg delivered years_since_LLFpre8-years_since_LLFpre2 years_since_LLFpost0-years_since_LLFpost7 ///
				pre_8_son-pre_2_son post_0_son-post_7_son ///
				p1_LLFpre8-p1_LLFpre2 p1_LLFpost0-p1_LLFpost7  ///
				p2_LLFpre8-p2_LLFpre2 p2_LLFpost0-p2_LLFpost7 ///
				p2_LLFpre8_son-p2_LLFpre2_son p2_LLFpost0_son-p2_LLFpost7_son ///
				p3_LLFpre8-p3_LLFpre2 p3_LLFpost0-p3_LLFpost7 ///
				p3_LLFpre8_son-p3_LLFpre2_son p3_LLFpost0_son-p3_LLFpost7_son ///
				p_1 p_2 p_3 hasson  p2_son p3_son $controls1 $controls2 $controls3 i.FE_t i.FE_prov ///
				if years_since_LLF<8 & rural==1, cluster(prov)
		}
		
		*****************************************************
		* Rural Only
		*****************************************************
		else if `spec' == 7 {
			xi: reg delivered years_since_LLFpre8-years_since_LLFpre2 years_since_LLFpost0-years_since_LLFpost7 ///
				pre_8_son-pre_2_son post_0_son-post_7_son ///
				p1_LLFpre8-p1_LLFpre2 p1_LLFpost0-p1_LLFpost7  ///
				p2_LLFpre8-p2_LLFpre2 p2_LLFpost0-p2_LLFpost7 ///
				p2_LLFpre8_son-p2_LLFpre2_son p2_LLFpost0_son-p2_LLFpost7_son ///
				p3_LLFpre8-p3_LLFpre2 p3_LLFpost0-p3_LLFpost7 ///
				p3_LLFpre8_son-p3_LLFpre2_son p3_LLFpost0_son-p3_LLFpost7_son ///
				p_1 p_2 p_3 hasson  p2_son p3_son culturalrev $controls1 $controls2 $controls3  i.FE_t i.FE_prov ///
				if year!=1967 & year !=1968 & years_since_LLF<8 , cluster(prov)
		}


	************************************************
	* (1.2) Preserve the Coefficients
	************************************************
	
		mat coef = e(b)
		mat temp = J(1,1, .) 
		
		local j = 1
		
		forval i = 2/8 {	
			local k = abs(`i' - 10)
			lincom p`j'_LLFpre`k' + pre_`k'_son 
			local est = r(estimate)
			
			mat temp1 = coef[1,"p`j'_LLFpre`k'"]
			mat temp2 = `est'
			mat temp = temp \ temp1 \ temp2
			}
			
		forval i = 9/16 {
			local k = abs(`i' - 9)
			lincom p`j'_LLFpost`k' + post_`k'_son 
			local est = r(estimate)
	
			mat temp1 = coef[1,"p`j'_LLFpost`k'"]
			mat temp2 = `est'
			mat temp = temp \ temp1 \ temp2
			}
		
		forval j = 2/3 {
		
		forval i = 2/8 {	
			local k = abs(`i' - 10)
			lincom p`j'_LLFpre`k' + pre_`k'_son + p`j'_LLFpre`k'_son
			local est = r(estimate)
	
			mat temp1 = coef[1,"p`j'_LLFpre`k'"]
			mat temp2 = `est'
			mat temp = temp \ temp1 \ temp2
			}
			
		forval i = 9/16 {
			local k = abs(`i' - 9)
			lincom p`j'_LLFpost`k' + post_`k'_son + p`j'_LLFpost`k'_son
			local est = r(estimate)
	
			mat temp1 = coef[1,"p`j'_LLFpost`k'"]
			mat temp2 = `est'
			mat temp = temp \ temp1 \ temp2
			}
		} 
					
		forval i = 2/8 {	
			local k = abs(`i' - 10)
	
			mat temp1 = coef[1,"years_since_LLFpre`k'"]
			mat temp = temp \ temp1
			}
			
		forval i = 9/16 {
			local k = abs(`i' - 9)
	
			mat temp1 = coef[1,"years_since_LLFpost`k'"]
			mat temp = temp \ temp1
			}
		
		mat robustness = robustness, temp
		}
	
	mat colnames robustness = "EventYear" "Parity" "Son" "Main" "cil" "ciu" "Spec1" "Spec2" "Spec3" "Spec4" "Spec5" "Spec6" "Spec7" 

		
	************************************************
	* (1.3) Plot Results
	************************************************
		
	preserve
		clear 
		svmat robustness, name(col)
			
		replace EventYear = EventYear - 9 
		replace EventYear = EventYear - 1 if EventYear <= -1
		
		forval i = 1/7 { 
			replace Spec`i'  = Spec`i'*100
		}
		
		replace cil = cil*100
		replace ciu  = ciu*100
		replace Main = Main*100
	
		drop if EventYear < -8 
		drop if EventYear > 8 

	gen yline = 0 
	replace Parity = Parity + 1
	

		*****************************
		* (1.3.2) Second Order Births
		*****************************
				
		replace EventYear = EventYear + .15 	if Son == 1
		replace EventYear = EventYear - .15 	if Son == 0
		
		twoway 	///
					(rcap cil ciu EventYear if Parity == 2 & Son==1, lcolor(blue))   ///
					(scatter Spec1 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec2 EventYear if Parity == 2 & Son==1, msymbol(circle_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec3 EventYear if Parity == 2 & Son==1, msymbol(diamond_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec4 EventYear if Parity == 2 & Son==1, msymbol(triangle_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec5 EventYear if Parity == 2 & Son==1, msymbol(x) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(rcap cil ciu EventYear if Parity == 2 & Son==0, lcolor(red))   ///
					(scatter Spec1 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec2 EventYear if Parity == 2 & Son==0, msymbol(circle_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec3 EventYear if Parity == 2 & Son==0, msymbol(diamond_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec4 EventYear if Parity == 2 & Son==0, msymbol(triangle_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec5 EventYear if Parity == 2 & Son==0, msymbol(x) mcolor(red) lcolor(red) cmissing(n)) ///
					(line yline EventYear if Parity == 2 & Son==., lcolor(black) cmissing(n) lpattern(solid)) ///
					(pci -24 -0.5 12 -0.5, lcolor(black)) ///
					,  ytitle("Percentage Point" "Change in Probability" ) ///
					xtitle("Years Since LLF Policy") yline(0, lcolor(black)) xline(-.5, lcolor(black)) xlab(-8(1)7)  ylab(-24(4)12) ///
					legend(off) title("Second Order Births")

			graph export `"$d_fig\Fig_A9a_PPR.jpg"', replace		
			
		twoway 	///
					(rcap cil ciu EventYear if Parity == 2 & Son==1, lcolor(blue))   ///
					(scatter Spec6 EventYear if Parity == 2 & Son==1, msymbol(circle_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec7 EventYear if Parity == 2 & Son==1, msymbol(diamond_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(rcap cil ciu EventYear if Parity == 2 & Son==0, lcolor(red))   ///
					(scatter Spec6 EventYear if Parity == 2 & Son==0, msymbol(circle_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec7 EventYear if Parity == 2 & Son==0, msymbol(diamond_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(line yline EventYear if Parity == 2 & Son==., lcolor(black) cmissing(n) lpattern(solid)) ///
					(pci -24 -0.5 12 -0.5, lcolor(black)) ///
					,  ytitle("Percentage Point" "Change in Probability" ) ///
					xtitle("Years Since LLF Policy") yline(0, lcolor(black)) xline(-.5, lcolor(black)) xlab(-8(1)7)  ylab(-24(4)12) ///
					legend(off) title("Second Order Births")

			graph export `"$d_fig\Fig_A9b_PPR.jpg"', replace	

		*****************************
		* (1.3.3) Third Order Births
		*****************************
	
		twoway 	///
					(rcap cil ciu EventYear if Parity == 3 & Son==1, lcolor(blue))   ///
					(scatter Spec1 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec2 EventYear if Parity == 3 & Son==1, msymbol(circle_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec3 EventYear if Parity == 3 & Son==1, msymbol(diamond_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec4 EventYear if Parity == 3 & Son==1, msymbol(triangle_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec5 EventYear if Parity == 3 & Son==1, msymbol(x) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(rcap cil ciu EventYear if Parity == 3 & Son==0, lcolor(red))   ///
					(scatter Spec1 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec2 EventYear if Parity == 3 & Son==0, msymbol(circle_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec3 EventYear if Parity == 3 & Son==0, msymbol(diamond_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec4 EventYear if Parity == 3 & Son==0, msymbol(triangle_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec5 EventYear if Parity == 3 & Son==0, msymbol(x) mcolor(red) lcolor(red) cmissing(n)) ///
					(line yline EventYear if Parity == 3 & Son==., lcolor(black) cmissing(n) lpattern(solid)) ///
					(pci -24 -0.5 12 -0.5, lcolor(black)) ///
					,  ytitle("Percentage Point" "Change in Probability" ) ///
					xtitle("Years Since LLF Policy") yline(0, lcolor(black)) xline(-.5, lcolor(black)) xlab(-8(1)7)  ylab(-24(4)12) ///
					legend(off) title("Third Order Births")

			graph export `"$d_fig\Fig_A10a_PPR.jpg"', replace
		
		twoway 	///
					(rcap cil ciu EventYear if Parity == 3 & Son==1, lcolor(blue))   ///
					(scatter Spec6 EventYear if Parity == 3 & Son==1, msymbol(circle_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec7 EventYear if Parity == 3 & Son==1, msymbol(diamond_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(rcap cil ciu EventYear if Parity == 3 & Son==0, lcolor(red))   ///
					(scatter Spec6 EventYear if Parity == 3 & Son==0, msymbol(circle_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec7 EventYear if Parity == 3 & Son==0, msymbol(diamond_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(line yline EventYear if Parity == 3 & Son==., lcolor(black) cmissing(n) lpattern(solid)) ///
					(pci -24 -0.5 12 -0.5, lcolor(black)) ///
					,  ytitle("Percentage Point" "Change in Probability" ) ///
					xtitle("Years Since LLF Policy") yline(0, lcolor(black)) xline(-.5, lcolor(black)) xlab(-8(1)7)  ylab(-24(4)12) ///
					legend(off) title("Third Order Births")

			graph export `"$d_fig\Fig_A10b_PPR.jpg"', replace


		*****************************
		* (1.3.4) Fourth Order Births
		*****************************
		
		twoway 	///
					(rcap cil ciu EventYear if Parity == 4 & Son==1, lcolor(blue))   ///
					(scatter Spec1 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec2 EventYear if Parity == 4 & Son==1, msymbol(circle_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec3 EventYear if Parity == 4 & Son==1, msymbol(diamond_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec4 EventYear if Parity == 4 & Son==1, msymbol(triangle_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec5 EventYear if Parity == 4 & Son==1, msymbol(x) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(rcap cil ciu EventYear if Parity == 4 & Son==0, lcolor(red))   ///
					(scatter Spec1 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec2 EventYear if Parity == 4 & Son==0, msymbol(circle_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec3 EventYear if Parity == 4 & Son==0, msymbol(diamond_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec4 EventYear if Parity == 4 & Son==0, msymbol(triangle_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec5 EventYear if Parity == 4 & Son==0, msymbol(x) mcolor(red) lcolor(red) cmissing(n)) ///
					(line yline EventYear if Parity == 4 & Son==., lcolor(black) cmissing(n) lpattern(solid)) ///
					(pci -24 -0.5 12 -0.5, lcolor(black)) ///
					,  ytitle("Percentage Point" "Change in Probability" ) ///
					xtitle("Years Since LLF Policy") yline(0, lcolor(black)) xline(-.5, lcolor(black)) xlab(-8(1)7)  ylab(-24(4)12) ///
					legend(off) title("Fourth Order Births")

			graph export `"$d_fig\Fig_A11a_PPR.png"', replace
		
		twoway 	///
					(rcap cil ciu EventYear if Parity == 4 & Son==1, lcolor(blue))   ///
					(scatter Spec6 EventYear if Parity == 4 & Son==1, msymbol(circle_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(scatter Spec7 EventYear if Parity == 4 & Son==1, msymbol(diamond_hollow) mcolor(blue) lcolor(blue) cmissing(n)) ///
					(rcap cil ciu EventYear if Parity == 4 & Son==0, lcolor(red))   ///
					(scatter Spec6 EventYear if Parity == 4 & Son==0, msymbol(circle_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(scatter Spec7 EventYear if Parity == 4 & Son==0, msymbol(diamond_hollow) mcolor(red) lcolor(red) cmissing(n)) ///
					(line yline EventYear if Parity == 4 & Son==., lcolor(black) cmissing(n) lpattern(solid)) ///
					(pci -24 -0.5 12 -0.5, lcolor(black)) ///
					,  ytitle("Percentage Point" "Change in Probability" ) ///
					xtitle("Years Since LLF Policy") yline(0, lcolor(black)) xline(-.5, lcolor(black)) xlab(-8(1)7)  ylab(-24(4)12) ///
					legend(off) title("Fourth Order Births")

			graph export `"$d_fig\Fig_A11b_PPR.jpg"', replace			
			
		restore	
				
*********************************	
* (2.0) "Leave One Out" Analysis
*********************************	
		
	*********************************	
	* (2.1) Regression
	*********************************	
	mat leaveoneout = results 
	tab prov, gen(prov_)
	local numprovs = r(r)
	
	forval province = 1/`numprovs' { 

	xi: reg delivered years_since_LLFpre8-years_since_LLFpre2 years_since_LLFpost0-years_since_LLFpost7 ///
				pre_8_son-pre_2_son post_0_son-post_7_son ///
				p1_LLFpre8-p1_LLFpre2 p1_LLFpost0-p1_LLFpost7  ///
				p2_LLFpre8-p2_LLFpre2 p2_LLFpost0-p2_LLFpost7 ///
				p2_LLFpre8_son-p2_LLFpre2_son p2_LLFpost0_son-p2_LLFpost7_son ///
				p3_LLFpre8-p3_LLFpre2 p3_LLFpost0-p3_LLFpost7 ///
				p3_LLFpre8_son-p3_LLFpre2_son p3_LLFpost0_son-p3_LLFpost7_son ///
				p_1 p_2 p_3 hasson  p2_son p3_son culturalrev  $controls1 $controls2 $controls3 i.FE_t i.FE_prov ///
				if year!=1967 & year !=1968 & years_since_LLF<8 & rural==1 & prov_`province'==0, cluster(prov)
	
	*********************************	
	* (2.2) Preserve the Coefficients
	*********************************	
		mat coef = e(b)
		mat temp = J(1,1, .) 
		
		local j = 1
		
		forval i = 2/8 {	
			local k = abs(`i' - 10)
			lincom p`j'_LLFpre`k' + pre_`k'_son 
			local est = r(estimate)
			
			mat temp1 = coef[1,"p`j'_LLFpre`k'"]
			mat temp2 = `est'
			mat temp = temp \ temp1 \ temp2
			}
			
		forval i = 9/16 {
			local k = abs(`i' - 9)
			lincom p`j'_LLFpost`k' + post_`k'_son 
			local est = r(estimate)
	
			mat temp1 = coef[1,"p`j'_LLFpost`k'"]
			mat temp2 = `est'
			mat temp = temp \ temp1 \ temp2
			}
		
		forval j = 2/3 {
		
		forval i = 2/8 {	
			local k = abs(`i' - 10)
			lincom p`j'_LLFpre`k' + pre_`k'_son + p`j'_LLFpre`k'_son
			local est = r(estimate)
	
			mat temp1 = coef[1,"p`j'_LLFpre`k'"]
			mat temp2 = `est'
			mat temp = temp \ temp1 \ temp2
			}
			
		forval i = 9/16 {
			local k = abs(`i' - 9)
			lincom p`j'_LLFpost`k' + post_`k'_son + p`j'_LLFpost`k'_son
			local est = r(estimate)
	
			mat temp1 = coef[1,"p`j'_LLFpost`k'"]
			mat temp2 = `est'
			mat temp = temp \ temp1 \ temp2
			}
		} 
					
		forval i = 2/8 {	
			local k = abs(`i' - 10)
	
			mat temp1 = coef[1,"years_since_LLFpre`k'"]
			mat temp = temp \ temp1
			}
			
		forval i = 9/16 {
			local k = abs(`i' - 9)
	
			mat temp1 = coef[1,"years_since_LLFpost`k'"]
			mat temp = temp \ temp1
			}
		
		mat leaveoneout = leaveoneout, temp
		}
	
	
	mat colnames leaveoneout = 	"EventYear" "Parity" "Son" "Main" "cil" "ciu" ///
								"prov1" "prov2" "prov3" "prov4" "prov5" "prov6" "prov7" "prov8" "prov9" "prov10" ///
								"prov11" "prov12" "prov13" "prov14" "prov15" "prov16" "prov17" "prov18" "prov19" "prov20" ///
								"prov21" "prov22" "prov23" "prov24" "prov25" "prov26" "prov27" "prov28"  

		
	*********************************	
	* (2.3) Plot Results
	*********************************	
	preserve
		clear 
		svmat leaveoneout, name(col)
			
		replace EventYear = EventYear - 9 
		replace EventYear = EventYear - 1 if EventYear <= -1
		
		forval i = 1/28 { 
			replace prov`i'  = prov`i'*100
		}
			
		replace cil = cil*100
		replace ciu  = ciu*100
		replace Main = Main*100
	
		drop if EventYear < -8 
		drop if EventYear > 8 

	gen yline = 0 
	replace Parity = Parity + 1
	

		****************************
		* (2.3.2) Second Order Births
		****************************		
		replace EventYear = EventYear + .15 	if Son == 1
		replace EventYear = EventYear - .15 	if Son == 0

		twoway 	///
				(scatter prov1 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov2 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov3 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov4 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov5 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov6 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov7 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov8 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov9 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov10 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov11 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov12 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov13 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov14 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov15 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov16 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov17 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov18 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov19 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov20 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///				
				(scatter prov21 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov22 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov23 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov24 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov25 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov26 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov27 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov28 EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(rcap cil ciu EventYear if Parity == 2 & Son==1,  lcolor(blue))   ///
				(scatter Main EventYear if Parity == 2 & Son==1, msymbol(circle) mcolor(blue) lcolor(blue) msize(small) cmissing(n)) ///
				(scatter prov1 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov2 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov3 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov4 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov5 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov6 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov7 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov8 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov9 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov10 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov11 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov12 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov13 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov14 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov15 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov16 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov17 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov18 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov19 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov20 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///				
				(scatter prov21 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov22 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov23 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov24 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov25 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov26 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov27 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov28 EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(line yline EventYear if Parity == 1 & EventYear>=-8, lcolor(black) cmissing(n) lpattern(solid)) ///
				(rcap cil ciu EventYear if Parity == 2 & Son==0,  lcolor(red))   ///
				(scatter Main EventYear if Parity == 2 & Son==0, msymbol(circle) mcolor(red) lcolor(red) msize(small) cmissing(n)) ///
				(pci -24 -0.5 12 -0.5, lcolor(black)) ///
				,  ytitle("Percentage Point" "Change in Probability" ) ///
				xtitle("Years Since LLF Policy") yline(0, lcolor(black)) xline(-.5, lcolor(black)) xlab(-8(1)7)  ylab(-24(4)12) ///
				legend(pos(6) rows(2) order(30 "Main Estimate, With Son(s)" 29 "95% Confidence Interval, With Son(s)" 61 "Main Estimate, Without Son(s)" 60 "95% Confidence Interval, Without Son(s)") ) title("Second Births")
			
		graph export `"$d_fig\Fig_A9c_PPR.jpg"', replace	
		
		****************************
		* (2.3.3) Third Order Births
		****************************		
	
		twoway 	///
				(scatter prov1 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov2 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov3 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov4 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov5 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov6 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov7 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov8 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov9 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov10 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov11 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov12 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov13 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov14 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov15 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov16 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov17 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov18 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov19 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov20 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///				
				(scatter prov21 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov22 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov23 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov24 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov25 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov26 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov27 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov28 EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(rcap cil ciu EventYear if Parity == 3 & Son==1,  lcolor(blue))   ///
				(scatter Main EventYear if Parity == 3 & Son==1, msymbol(circle) mcolor(blue) lcolor(blue) msize(small) cmissing(n)) ///
				(scatter prov1 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov2 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov3 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov4 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov5 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov6 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov7 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov8 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov9 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov10 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov11 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov12 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov13 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov14 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov15 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov16 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov17 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov18 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov19 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov20 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///				
				(scatter prov21 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov22 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov23 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov24 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov25 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov26 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov27 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov28 EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(line yline EventYear if Parity == 3 & EventYear>=-8, lcolor(black) cmissing(n) lpattern(solid)) ///
				(rcap cil ciu EventYear if Parity == 3 & Son==0,  lcolor(red))   ///
				(scatter Main EventYear if Parity == 3 & Son==0, msymbol(circle) mcolor(red) lcolor(red) msize(small) cmissing(n)) ///
				(pci -24 -0.5 12 -0.5, lcolor(black)) ///
				,  ytitle("Percentage Point" "Change in Probability" ) ///
				xtitle("Years Since LLF Policy") yline(0, lcolor(black)) xline(-.5, lcolor(black)) xlab(-8(1)7)  ylab(-24(4)12) ///
				legend(pos(6) rows(2) order(30 "Main Estimate, With Son(s)" 29 "95% Confidence Interval, With Son(s)" 61 "Main Estimate, Without Son(s)" 60 "95% Confidence Interval, Without Son(s)") ) title("Third Births")
			
		graph export `"$d_fig\Fig_A10c_PPR.jpg"', replace	

		****************************
		* (2.3.4) Fourth Order Births
		****************************		
		
		twoway 	///
				(scatter prov1 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov2 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov3 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov4 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov5 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov6 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov7 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov8 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov9 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov10 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov11 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov12 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov13 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov14 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov15 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov16 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov17 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov18 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov19 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov20 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///				
				(scatter prov21 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov22 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov23 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov24 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov25 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov26 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov27 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov28 EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///		
				(rcap cil ciu EventYear if Parity == 4 & Son==1,  lcolor(blue))   ///
				(scatter Main EventYear if Parity == 4 & Son==1, msymbol(circle) mcolor(blue) lcolor(blue) msize(small) cmissing(n)) ///
				(scatter prov1 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov2 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov3 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov4 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov5 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov6 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov7 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov8 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov9 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov10 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov11 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov12 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov13 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov14 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov15 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov16 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov17 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov18 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov19 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov20 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///				
				(scatter prov21 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov22 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov23 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov24 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov25 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov26 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov27 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(scatter prov28 EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(black) lcolor(black) cmissing(n)) ///
				(line yline EventYear if Parity == 4 & EventYear>=-8, lcolor(black) cmissing(n) lpattern(solid)) ///
				(rcap cil ciu EventYear if Parity == 4 & Son==0,  lcolor(red))   ///
				(scatter Main EventYear if Parity == 4 & Son==0, msymbol(circle) mcolor(red) lcolor(red) msize(small) cmissing(n)) ///
				(pci -24 -0.5 12 -0.5, lcolor(black)) ///
				,  ytitle("Percentage Point" "Change in Probability" ) ///
				xtitle("Years Since LLF Policy") yline(0, lcolor(black)) xline(-.5, lcolor(black)) xlab(-8(1)7)  ylab(-24(4)12) ///
				legend(pos(6) rows(2) order(30 "Main Estimate, With Son(s)" 29 "95% Confidence Interval, With Son(s)" 61 "Main Estimate, Without Son(s)" 60 "95% Confidence Interval, Without Son(s)") ) title("Fourth Births")
			
		graph export `"$d_fig\Fig_A11c_PPR.jpg"', replace	
		
		
	restore	
	